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I. INTRODUCTION 

Due to fast experimental progress [IH1] , it becomes in- 
creasingly important to model quantum optical systems 
with spontaneous photon emission more and more accu- 
rately [5H7]. Motivated by this, we revisit the derivation 
of the master equation of an atomic two-level system, like 
a single quantum dot, a single colour centre in diamond, 
or a well-localised atom or ion, within a free radiation 
field without employing the rotating wave approxima- 
tion (RWA) and certain Markovian approximations. The 
starting point of our derivation is, as usual, a quantum 
electrodynamic Hamiltonian of the form 

H = H A +H F + H AF , (1) 

where Ha, Hp and Haf are the atom, the field, and 
the atom-field interaction components. This Hamilto- 
nian acts on a Hilbert space % — Ha ® %f whose sub- 
spaces Ha and Hf are the eigenspaces of Ha and Hp, 
respectively. The purpose of the master equation is to 
summarise the effective time evolution of the atomic sys- 
tem induced by this Hamiltonian and by the coupling of 
the field to a photon-absorbing environment. 

But there is a problem: There are in principle in- 
finitely many unitarily equivalent representations H' of 
the Hamiltonian H in Eq. ([I]) . Each of these Hamiltoni- 
ans relates to H via a unitary transformation R, 

H' = RHRK (2) 

Now, suppose R does not commute with Ha and Hp. 
Then, although physically equivalent, H and W corre- 
spond to different decompositions of the Hilbert space % 
into subspaces Ha and Hf ■ Both refer in general to phys- 
ically different components as representing the "system" 
and as representing the "field" [8]. The Hamiltonians 
H and H' therefore result in general in different mas- 
ter equations. This problem only vanishes when certain 
Markovian approximations are employed which reduce 
every master equation to the standard Born-Markov form 
- independent of the chosen representation of H . When 



going beyond these approximations, we are faced with 
the question as to which representation of H should be 
employed? 

That there are many different representations of H is 
directly related to the fact that there are many possible 
choices of gauge of the electromagnetic field. The most 
prominent form of H is the so-called minimal coupling 
Hamiltonian which is obtained after applying the canon- 
ical quantisation procedure to the classical Lagrangian 
in Coulomb gauge [3]. Starting from this Hamiltonian, a 
more general form of the atom-field Hamiltonian, includ- 
ing the multipolar Hamiltonian, is obtained through the 
application of the Power-Zienau-Woolley (PZW) trans- 
formation [TtilfT2"] . Further representations of H are ob- 
tained through a generalised PZW transformation which 
depends on a set of parameters a^. As pointed out by 
Drummond [5] , a particular choice of au 's can be made to 
eliminate the counter-rotating terms in the atom-field in- 
teraction Haf > when assuming that the size of the atomic 
system is much smaller than a typical optical wave length 
(electric dipole approximation) . In the following we refer 
to the resulting Hamiltonian as the rotating wave Hamil- 
tonian, since almost the same Hamiltonian, ie. a Hamil- 
tonian of the same form but with different coupling con- 
stants, arises when applying the RWA. 

This paper considers the vast set of unitarily equiv- 
alent atom-field Hamiltonians generated by the above 
mentioned generalisation of the PZW transformation [8] 
and derives a master equation for every one of them. It 
is then shown that most of these equations predict enor- 
mous stationary state photon emission rates I ss , even in 
the absence of external driving, when we consider the pa- 
rameters of recent experiments with single quantum dots 
[3] and single colour centers in diamond [I]. The only 
master equation which yields I ss = is the one associ- 
ated with the rotating wave Hamiltonian. This means, 
this Hamiltonian (and not the minimal coupling or the 
multipolar Hamiltonians) constitutes the physically most 
relevant representation of the atom-field Hamiltonian in 
Eq. 0. As we shall see below, the corresponding master 
equation is essentially the same as the standard Born- 
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Markov master equation. However, it is no longer the re- 
sult of approximations but the result of a specific choice 
of R which identifies the bare states of the atomic system 
and its free radiation field correctly. 

Many different ways of deriving master equations, ie. of 
deriving an effective equation for the time evolution of 
the atomic system, can be found in the literature. Due 
to a great variety of microscopic models, different au- 
thors usually apply different approximations. For exam- 
ple, many derivations of master equations rely on the 
so-called Born approximation [13-18]. This approxima- 
tion assumes that the state of a system and its bath can 
be written as a product state of the form pa <8> Pf at 
any possible initial time t, where pp denotes the thermal 
equilibrium of the bath. In general it is emphasized that 
the Born approximation applies well only in certain sit- 
uations, like the weak coupling regime, where memory 
effects remain negligible [HI [20] . 

Moreover, many authors emphasize that there is a 
strong overlap between the validity range of the Born 
approximation and the validity range of certain Marko- 
vian approximations. A system is Markovian, ie. with- 
out a memory, when correlations between system and 
environment exist only on time scales which are much 
shorter than the time scale on which the system evolves. 
It has therefore been argued that it is not possible to 
apply the Born approximation without using in addition 
certain Markovian approximations |18j . When applying 
both approximations to atom-field systems, we obtain the 
standard Born-Markov master equation [2T] . All possible 
representations of H result in the same master equation, 
independent of whether or not the RWA is applied to the 
Hamiltonian. In other words, the only way of obtain- 
ing a master equation which is more precise than already 
existing master equations is to allow for non-Markovian 
descriptions (cf. eg. Refs. [22H24] 1 . However, although 
these papers contain interesting results for a variety of 
quantum systems, they do not apply to the system which 
we consider here. 

The reason for this is that the above observation only 
applies when system and bath form a closed quantum 
system with only unitary dynamics, as it is widely as- 
sumed in the literature fBHTSI |2"5] . Here we adopt a 
different point of view. Following the ideas of Hegerfeldt 
and Wilser [13 [27] and others HBHH], we take the in- 
teraction of the free radiation field with an external en- 
vironment, like a detector or the walls of the laboratory, 
explicitly into account. Following the ideas of Zurek and 
Paz on decoherence and einselection (see eg. Ref. [3T]) 
and applying them to the free radiation field and its ex- 
ternal environment, we assume that this environment re- 
sets the free radiation field on a coarse grained time scale 
At onto its vacuum state |0). This means, the Born ap- 
proximation which assumes that the density matrix of 
the atom-field system can be written as 

P(t) = p A (t)®|0><0|, (3) 

is no longer an approximation but the result of external 



environmental couplings. Hence it is well justified to de- 
rive an effective master equation for the atomic system 
without employing the RWA and certain Markovian ap- 
proximations, but still employing Eq. ([3]). As a result, 
we obtain a master equation whose constants depend on 
the typical environmental response time At and which is 
automatically Markovian. This is as expected, since the 
free radiation field does not constitute a memory for the 
states of the atomic system. 

The above described approach of deriving master equa- 
tions for atomic systems with spontaneous photon emis- 
sion is in very good agreement with actual quantum 
optics experiments |32H35j . For example, it has been 
shown experimentally that certain three-level atoms ex- 
hibit macroscopic quantum jumps |32H34| . ie. they emit 
a completely random telegraph signal with long intervals 
of intense photon emissions interrupted by periods of the 
complete absence of photons [361 137j . In the experiment 
by Eichmann et al. [35] two laser driven two-level atoms 
emit photons onto a far-away screen, where they form 
an interference patterns which is typical for double slit 
experiments |38) . In both experiments, the free radiation 
field constantly leaks information about the atomic sys- 
tem into an external environment. When an atom emits 
a real photon, the photon flies away and is absorbed after 
a relatively short time At. It does not return to interact 
again with its source. 

As summarised in Eq. ([3]), environmental couplings 
project the atom-field system onto a product state with 
the free radiation field in its vacuum state. Notice how- 
ever that for most representations of the atom-field sys- 
tem, the ground state of the Hamiltonian H in Eq. ([I]) 
is an entangled state. This means, in most derivations 
of master equations, the interaction with the environ- 
ment resets the atom-field system constantly onto a state 
which is energetically higher than its ground state. If this 
applies, the photon-absorbing environment constantly 
pumps energy into the system which then manifests itself 
as a non-zero stationary state photon emission rate even 
in the absence of external driving [5S] 2D] . Such a con- 
centration of energy might play a crucial role in sonolu- 
minescence experiments [41j . However, it need not play a 
role in atom-field systems, since the ground state of the 
rotating wave Hamiltonian is a product state with the 
free radiation field in its vacuum state. 

Some authors argue that the description of a physi- 
cal system like an atom in a free radiation field should 
be gauge-independent (see eg. Ref. [12] and references 
therein). This is no longer the case when we select for 
example the master equation associated with the rotat- 
ing wave Hamiltonian. However, notice that we consider 
the atom-field system as an open quantum system with a 
photon- absorbing environment. This environment makes 
a concrete choice as to which component of the atom-field 
system represents the "field" and which one represents 
the "atom." A similar choice should be made when de- 
riving an atomic master equation. 

There are five sections in this paper. In Section [n] 
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we introduce a large class of unitarily equivalent repre- 
sentations of the atom-field Hamiltonian. In Section IIII 



we obtain a large set of non-equivalent master equations. 
In Section |TV| we calculate the corresponding stationary 
state photon emission rate. The results of this calculation 
help us to select what we believe to be the most accurate 
master equation. Finally, we summarise our findings in 
Section M 



II. UNITARILY EQUIVALENT HAMILTONIANS 

In this section we focus our attention on a large set 
of unitarily equivalent atom-field Hamiltonians which we 
obtain by applying a wide range of unitary transforma- 
tions R to the minimal coupling Hamiltonian. For ex- 
ample, the Power- Zienau-Woolley (PZW) transformation 
transforms this Hamiltonian into the so-called multipolar 
Hamiltonian [rUlfi"2"] . In general, we obtain a Hamilto- 
nian which is a mixture of the minimal coupling and the 
multipolar Hamiltonians [33] . The most notable of these 
Hamiltonians is the one that does not contain any counter 
rotating terms in the atom-field interaction 8] . Since this 
Hamiltonian is of the same form as the Hamiltonian H of 
the atom-field system obtained when applying the RWA, 
we refer to it in the following as the rotating wave Hamil- 
tonian. This Hamiltonian is now no longer the result of 
approximations but the result of a certain unitary trans- 
formation R. 



A. The minimal coupling Hamiltonian 

To begin with we consider an electron of charge e and 
mass m which is trapped inside a potential V(r) and that 
couples to a continuum of quantised electromagnetic field 
modes. The canonical operators of the radiation field 
are the transverse vector potential A and its conjugate 
momentum II. These satisfy the transverse canonical 
commutation relation 



[A^H^x')] = iS^(x-x') 



(4) 



where the indices i and j specify different spatial com- 
ponents. The canonical operators for the electron are 
its position r and its conjugate momentum p with the 
canonical commutation relation 

[n,pj] = ihSij . (5) 

Using this notation, the minimal coupling Hamiltonian 
-ffmin can be written as 
1 
2m 
1 



■ eA(O)] 
1 



d J x 



(4) 



' + V(t) 
B(x)| 2 - 



1 



n(x) 



(6) 



The integration in this equation is over the whole space 
occupied by the atom-field system and B = V x A de- 
notes the magnetic field. The first term on the right hand 



side of Eq. ^ has already been simplified with the help 
of the electric dipole approximation which replaces the 
position operator r of the electron in the vector poten- 
tial A with the position of the centre of mass of the 
particle. 

One can easily check that two operators A and II de- 
fined as 



A « = E 



kA 



\2eoUJkL 3 



( heoUk \ 



H.c. 



n(x) = -iE(^r e ^kAC ik - x + H.c. (7) 



kA 



V 2L3 J 



obey the canonical commutator relation in Eq. Q , if the 
a kA arid <zt» are creation and annihilation operators with 



i kA 

the bosonic commutation relation 



[ftkAiffljL/A'] = <^kk' £> 



AA' 



(8) 



and with an appropriate choice of boundary conditions. 
Here L 3 denotes the quantisation volume of the free ra- 
diation field. Each field mode (k, A) is characterised by 
a polarisation A, a wavevector k, and orthogonal unit 
polarisation vectors e kA with e k i, e k 2, and k all being 
pairwise orthogonal. The corresponding magnetic field 
B is in addition given by 



B(x) = i 



-H.c, 



ik-x 



(k x e k A)a k Ae 1 ' 



(9) 



where k = k/fc is a normalised vector. Eqs. (|7|)-(|9|) can 
now be used to bring the minimal coupling Hamiltonian 
into a more compact form, 



tfmi„ = — [p + eA(0)] 2 + V^(r) + ^ 7^,. a^kA 

kA 



(10) 

where only an overall constant with no physical conse- 
quences has been neglected. 



B. Alternative representations of the atom-field 
Hamiltonian 

As suggested by Drummond [5], we now consider the 
unitary transformation R{ ak } in the electric dipole ap- 
proximation which is defined as 



R {ak} = exp ( -A W} (0) -r 



(11) 



The field operator A{ afc }(0) in this equation is given by 

i 

A W}(°) = E ( o, . .. r J a fe e kA a kA + H.c. (12) 



kA 



2e Lu k L 3 
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with the cafe's being real and dimensionless. In the case, 
where all the otk in Eq. (11 1 equal one, the above trans- 
formation R{a k } reduces to the well-known PZW trans- 
formation in the electric dipole approximation [SJ 03] . 

We now apply R{a k } in Eq. (Illl) to the minimal cou- 
pling Hamiltonian in Eq. ( [To] ) . The details of this cal- 
culation which yields an {a^j-dependent Hamiltonian of 
the same form as Eq. ([I]) but with its components given 

by m 



H F = 



P 

2m 
2m 



2e L3 



A(0) - A {Qfc} (0)| 2 + hoj k 4 A a kA 



k,\ 



p • (A(0) - A {Qfc} (0)) - - n {Qfe} (0) • r 



(13) 



can be found in App. IaI The field operator Ur ak y(0) in 
this equation is given ay 

i 

n W}(°) = - i E( fi 2Z?) 2 «*e kA o kA +H.c. (14) 

The transformation does not leave the different 

components of the Hilbert spaces Ha and Hf of atom 
and field invariant. This is because the canonical mo- 
menta p and II take on different physical meanings for 
each one of the possible representations of the atom-field 
Hamiltonian H . Each representation refers to physically 
different components of the total system as representing 
the "atom" and as representing the "field." 



The bare Hamiltonians Ha and Hp in Eq. ( 13 1 contain 
self-energy contributions. These are divergent without 
the introduction of an upper frequency cut-off for the 
modes of the free radiation field. In order to show that 
these divergences do not contribute to our derivation of 
master equations, we now consider an atomic two-level 
system with ground state |1), excited state |2), and an 
energy separation Hujq. Moreover, we introduce the Pauli 
operators 



|2)(1|, ( r- = |l)(2|, 
cxs = | (|2><2! - ]1><1|) 

with the spin-like commutation relations 

[a + ,a-]=2a 3 , [a ± , a 3 ] = ■ 



Using the commutation relation 

[H A ,n] -- 



Pi 



(15) 



(16) 



(17) 



and multiplying the atomic operators r and p on both 
sides with the identity |1)(1| + \2)(2\ = 1, wc find that 
the components Ti and Pi may be written as 



n = di <r + + H.c. , 
Pi = imujQ di a + + H.c. 



(18) 



where the di are the components of the (in general com- 
plex) atomic dipole moment 



d= (2|r|l) 



(19) 



For a two-level atom, the atomic self-energy term merely 
shifts the zero-point energy, because it is diagonal in the 
bare atomic basis. Neglecting it is equivalent to neglect- 
ing a constant in the Hamiltonian. The field self-energy 
term does not contribute in the second order perturba- 
tion theory calculations in Section |IIID| It can therefore 
be neglected without affecting the results of this paper. 
Taking this into account, we obtain the general atom-field 
Hamiltonian 

# = E ^3kA c + (u^ al x + u% a kA ) + H.c. 

kA 

+fkv cr 3 + tiujk Q kA akA (20) 

kA 

with the atom-field coupling constant g k A and the (real) 
coefficients u k defined as 



.9kA 



(1-afc) 



^0 



1/2 



UJk_ 
UJ 



1/2 



(21) 



When setting all equal to zero, the minimal coupling 
Hamiltonian in the electric dipole approximation reduces 
to the minimal coupling form 



H„ 



E 

kA 



%kA 

+Huj a 3 



\ 1/2 



OkA 
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k a kA a k A 



(22) 



kA 



which contains a linear coupling between its respective 
atomic system and the modes (k, A) of its respective free 
radiation field. Let us now have a closer look at two other 
examples of the general Hamiltonian. 



C. The multipolar Hamiltonian 

Another Hamiltonian which is often used in the liter- 
ature is the multipolar Hamiltonian. It arises from the 
minimal coupling Hamiltonian after applying the PZW 
transformation. This means, the multipolar Hamilto- 
nian for an atomic two-level system and in the electric 
dipole approximation corresponds to the Hamiltonian H 
in Eq. ( 20 1 with all a% identical to one. It equals 



E^kA 



kA 



1/2 



+ C ) UkA 



4 



+HUJ Q <7 3 + ^ «kA a kA • 



(23) 



kA 
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Like the minimal coupling Hamiltonian, this represen- 
tation implies a linear coupling between its respective 
atomic system and the modes (k, A) of its respective free 
radiation field. 



D. The rotating wave Hamiltonian 

The physical motivation for introducing the unitary 
transformation Ri ak } m Eq. (11) is that this transfor- 
mation allows us to remove the counter rotating terms of 



the form u^a^x and er+aj^ from the atom-field Hamilto- 
nian in electric dipole approximation 8 . Indeed, when 
choosing the coefficients a k such that 



Ok 



W 



LOq + LO k 



(24) 



the general Hamiltonian H in Eq. ( 20 I simplifies to 



%kA — CT OkA + H.C 



kA 



LOq + UJ k 



kA 



(25) 



In the following, we refer to this Hamiltonian as the ro- 
tating wave Hamiltonian, since almost the same Hamil- 
tonian is obtained when applying the RWA. The only 
difference between these Hamiltonians is a difference in 
the atom-field coupling constants. 



this Hamiltonian into the interaction picture with respect 
to the free Hamiltonian 



Hq = HuJq <T 3 + ^ akA a kA 



(26) 



kA 



Using Eq. (20), we find that the general atom-field inter- 



action Hamiltonian Hi equals 



kA 

+H.c. 



-iu)i.t \ 
) 

(27) 



Within this interaction picture, the atomic density ma- 
trix px{t) in the Schrodinger picture becomes 

PA i(t) = U^(t,0)p A (t)U Q (t,0). (28) 
Defining the superoperator C such that 

p A (t) = C( PA (t)) (29) 



and taking the time derivative of Eq. ( 28 1 , we find that 
p A i{t) = iojq [c3,PAi] 

+C/ o t (i,0)£(c/ o (t,0)pAiC/ o t ( < :0)) U (t,0). (30) 

In the following subsections, we calculate this time 
derivative using second order perturbation theory and 
then use Eq. ( 30 ) to find the superoperator C in the 



Schrodinger picture. 



III. NON-EQUIVALENT MASTER EQUATIONS 

In this section we derive the the master equation for the 
effective time evolution of an atomic system with spon- 
taneous photon emission without applying the Born and 
certain Markovian approximations. Our only approxi- 
mations are the use of second order perturbation theory 
within an appropriately chosen interaction picture and 
the resetting of the free radiation field on a coarse grained 
time scale At into its vacuum state. The first assump- 
tion is well justified as long as the environmental response 
time At is short compared to the characteristic time scale 
of the effective atomic evolution. The second approxima- 
tion is justified in the presence of a photon-absorbing en- 
vironment which monitors the free radiation field [27H30) . 
Independent of whether or not a photon has been found, 
decoherence destroys any coherences between the radia- 
tion field and its surroundings and transfers it into an 
environmentally preferred (pointer or einselected) state 
- the vacuum state |0) [3"T] . 

A. Interaction picture 

The starting point of our derivation is, as usual, the 
Hamiltonian H in Eq. (TlJ) . In this subsection, we transfer 



B. Photon-absorbing environment 

Following ideas of Hegerfeldt and Wilser (5^1 [STJ and 
others [28, 29J, we do not assume that the atomic sys- 
tem and the surrounding free radiation field are a closed 
quantum system with purely unitary dynamics. Instead, 
we take the coupling of the free radiation field to an addi- 
tional external environment, like a detector or the walls 
of the laboratory, explicitly into account. This external 
environment thermalises very rapidly and is in general in 
an equilibrium state. We denote this equilibrium state 
in the following by p™ v . 

Suppose, the interaction between the atomic system 
and the free radiation field had already created a small 
photon population in the field modes (k, A). If the pho- 
tons in the free radiation field are real photons, they fly 
away and reach the external environment, ie. a detec- 
tor or the walls of the laboratory, after a relatively short 
time. There they interact with a large collection of atoms 
such that 

PF ®Ptr — ► io)(oi®p° nv ( PF ), (3i) 

where p-p is the initial density matrix of the free radia- 
tion field and |0) is its vacuum state with a^x 1 0) = 
for all k and A. In other words, any initial excitation 
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in the free radiation field vanishes very quickly as a re- 
sult of the interaction between the radiation field and its 
external environment. Eq. (31) also assumes that there 



are no coherences between the free radiation field and 
its environment. These have already been destroyed by 
decoherence [3T] . 

Moreover, the environment itself thermalises in general 
very rapidly Taking this into account, the effect of the 
photon-absorbing environment can be summarised as 



Pf 



|o><o|®pS 



(32) 



if the free radiation field does not become re-populated in 
the process. This applies when the vacuum state |0) is an 
einselected state of the free radiation field. Einselectcd 
(ie. environmentally-selected) states are states which re- 
main stable in spite of the environment [3"T] . Here we 
assume that the free radiation field possesses only a sin- 
gle einselected state. This is justified on the grounds 
that emitted photons fly away from the atom and are 
absorbed. In the absence of an atomic emitter but in the 
presence of a photon-absorbing environment, the vacuum 
state |0) is the only state of the free radiation field which 
does not evolve in time. 



C. Derivation of master equations 

The above considerations allow us to write the density 
matrix p(t) of the atom-field system at any time t as 



Pi{t) = PAi(t) 



(33) 



where pAi(t) denotes an atomic density matrix with 
respect to the above introduced interaction picture. 
Eq. ( |33[ ) is known as the Born approximation but, as 
we have seen above, the assumption of an uncorrelated 
atom-field state becomes exact in the presence of a 
photon- absorbing environment. 

The second step in our derivation of a master equation 
for the atomic system is the calculation of the time evo- 
lution of the state in Eq. (33) over a time interval At. 



Within this time, p\{t) in Eq. (33) evolves into 



At) 

U I (t + At,t)p M (t)® \0)(0\ Ul(t + At,*). (34) 



This atom-field density matrix corresponds in general to 
an entangled state with population in most modes (k, A) 
of the free radiation field. Combining Eqs. (31), (32), 



and (34), we find that the density matrix of atom, field, 



and environment evolves during any subsequent photon 
absorption according to 



Tr P 



Ui (t + At, t) p Al (t) ® |0) <0| U} (* + At, t) 



m^\®pt 



(35) 



where Trp denotes the trace over all the modes of the 
free radiation field. This equation shows that the density 
matrix of atom and field remains uncorrelated on the 
coarse grained time scale given by A*. 



An alternative way of deriving Eq. ( 35 ) is to assume an 



environment which performs rapidly repeated, photon- 
absorbing measurements on a coarse grained time scale 
A* [1H1HI]- Calculations based on the RWA [38], which 
reproduce the observations of recent quantum optics ex- 
periments like the one by Eichmann et al. [35], show 
that it is sufficient to assume an environment which re- 
solves the direction of each emitted photon but not their 
frequency in order to obtain Eq. ( 35 ) . Our derivation 



of master equations is hence consistent with the work 
by Hegerfeldt and Wilser [27] who derived master 
equations while assuming photon-absorbing broadband 
measurements. Like einselection, these environment- 
induced measurements destroy any correlations between 
the atomic system and the free radiation field. The trace 
over the free radiation field in Eq. (35) assures that the 



density matrix on the right hand side is always nor- 
malised. Eq. ( [35] ) also takes into account that the interac- 
tion between radiation field and its external environment 
occurs on a very short time scale without non-local effects 
on the atomic state. 

Here we are only interested in the time evolution of 
PAi(t)- Eq. (35) shows that this density matrix evolves 
such that 

PM (t + At) = p° A1 (t + At) + p^it + At) , (36) 



if we define 
p° AI (t + At) 

P^(* + At) 



(0|C/i(t + Ai,*)|0) pAi{t) 

x (0\ U}{t + At, *)|0), 

OO 

]T 5> kA | CM* + At,*) |0) paiW 



71=1 kA 



x <0| U}(t + At, *)|n kA ). (37) 



These atomic density matrices, respectively, describe the 
subensemble of atoms without and the subensemble of 
atoms with photon emission in (t, t + At). 

The purpose of an atomic master equation is to estab- 
lish a direct connection between p A \(t + At) in Eq. (36) 



and PAi(t) hi Eq. (33) without having to calculate the 
trace over the free radiation field after each time step 
A*. To create such a link, we calculate the difference 
quotient 



A(paiW) 



PAl(t + At) - P Al(t) 

At 



(38) 



In the limit A* — > 0, this expression becomes a time 
derivative, giving the relation 



PA\{t) = Ci{p A i(t)) 



(39) 



However, since the superoperator C\ is obtained via a 
partial trace, by applying this limit one enters a regime 
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in which the dynamics of the atom-field system becomes 
restricted onto the zero photon subspace [27]. We there- 
fore avoid this limit and only assume that the time in- 
crement At is sufficiently small compared to the effective 
time scale over which the atomic system evolves. Eq. ( 38 ) 



now constitutes a Markovian approximation of the con- 
tinuous dynamics given by Eq. ( 39 1 . 



D. Second order perturbation theory 

The next step in our derivation of an atomic master 
equation, ie. of the superoperator C in Eq. (291, is the 



explicit calculation of C\ in Eq. ( 39 ) . Since the time evo- 



lution of the atomic system is homogeneous in time, C\ 
is independent of t as long as no time-dependent inter- 
actions, like laser fields, are applied. However, one can 
easily check that C\ depends explicitly on At. In the fol- 
lowing, we treat At as an external parameter which can 
be absorbed into the coefficients of the relevant master 
equation. 



1. The subensemble without photon detection 

As mentioned already above, here we are especially 
interested in the case where At is small compared to 
the time scale on which the atomic density matrix pAi(t) 
evolves. This assumption allows us to calculate p%j(t + 
At) in Eq. (37) using second order perturbation theory. 



Doing so we find that 
(0\Uj(t + At,t)\0) 



1 



1 



t+At 



At' / dt"(0|ffi(i')W)|0). (40) 



Using Eq. (27), (a ) — 0, and introducing the time 
independent constants A' ± as 



1 r At r 



2 1( ± 2 p Ti("o±"fc)(*'-*") 



the above expression simplifies to 



(41) 



<0|l7i(* + At,t)|0) 

= 1 - (A'_a+a- +A' + ct-(t + ) At. (42) 



The constants A'± in Eq. (41) are in general complex 
numbers. Since their real and their imaginary parts have 
different physical meanings, it is convenient to consider 
them separately and to define 



A ± = 2Re(A' ± ). 
Awi, 2 = Im(A' ± ) 



(43) 



such that A'± = \A± + i Acj^ 2 . Substituting Eq. (40) 
into Eq. (37) and using this notation, we find that 



p° AI {t + At) = p AI (t)+iJ2^[pAi(t),\i)(i\}^t 

i=l 

-i (A-a+a- + A + <j-<r + )p AI (t) At 

--pAi{t) {A-a+cT- + A + <j-ct + ) At (44) 

up to first order in At. This means, the imaginary parts 
of A' + and A'_ shift respectively the frequencies of the 
atomic states |1) and |2), while the real parts of these 
constants are decay rates. 



2. The subensemble with photon detection 

Evaluating U\{t + At, t) again using second order per- 
turbation theory, the density matrix p A i(t + At) in 
Eq. (37) for the subensemble with photon emission in 



(t + At, t) becomes 

P> A1 (t + At) = / dt> 

l„ \ J t J t 



t+At 



At" 



(IkAl^l(t')IO) PAi(t) (0|fli(O|lkA> • (45) 

Substituting Eq. (27) into this equation, p'p il (t + At) sim- 
plifies to 



pl^t + At) = TZj(p A i(t))At, 

having defined the reset superoperator IZi as 

^l(PAl) = A + <T + p A i(J~ + A- er~pAic + 
+B + ct + pai(t + + B- a~ pA\o~ 



with 



A± = J2^k\\ 2 u^ Atsinc 2 

kA 

where sine a; = sinx/x, and with 

rt+At r t+At 



{uj ±ujk)At 



(46) 



(47) 



(48) 



B, = 



Jt Jt kA 



5k A «fc «fc 



i nt-\~At pt-\-L\Z 



X gi^o(*'+t") gi^fe(t'-t") 
t+At 



kA 

xe -iu (t'+t") e iw fc (t'-t") _ 



Using the specific properties of trigonometric functions, 
one can show that the A± in this equation are twice the 
real parts of the constants A'± in Eq. (41 ) and hence the 
same as the A± which we introduced in Section HID 
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FIG. 1: Logarithmic plot of the transition rate A- in Eq. (B3 1 FIG. 2: Logarithmic plot of the transition rates A+ in 



for the minimal coupling, the multipolar and the rotating 
wave Hamiltonian as a function of the upper cut-off frequency 
Umaj,, while LOmin = and o;oAi = 10 4 . 



3. The superoperator L\ 



Substituting Eqs. (44 1, (46 1, and |47| into Eq. (38), we 




find that the superoperator C\ in Eq. (|38|) equals 

2 



A(pai) 



i^Z Auj i [/°ai 5 

=i 

-A + (2a + p A ia~ - a~a + p A i - p A ia~a + ) 



i=i 
1 

1 



-A- (2a paict - v o PAi - PAW v 



2 

-B + a + p A icr + + B_ a~p A ia~ 



(50) 



The only approximation made in the derivation of this 
equation is the use of second order perturbation theory. 



E. General master equation 

To determine the superoperator C in the Schrddinger 
picture, we use Eq. ( 30 1 which implies 



£(pa) = -iwo [03 , pa] 

+U Q (t,0)& (ul(t,0)p A U o {t 1 0)) £/ o f (i,0). (51) 



Combining this equation with Eq. ( 50 1 and redefining the 



atomic transition frequency uiq such that the first term 
in Eq. (50) can be absorbed into the free energy of the 



atom, we are finally able to calculate the superoperator 



Eq. ( B3 1 for the minimal coupling and the multipolar Hamil- 
tonian as a function of the upper cut-off frequency o> max , while 
= and u) At = 10 4 . These are in very good agreement 



with the analytical results in Eqs. (B6HB8) 



£. As a result we obtain the general master equation 



Pa 



-uj [cr 3 , p A ] 

+ \ A + ( 2<j+ Pav~ - <r~a + p A - Pao-~ct h 



1 



A_ (2a p A a + — a + a p A — p A a^ a 



2 

+B a + pA<J + + B* a~ p A a~ 



(52) 



where wo is a shifted atomic frequency and where B is 
defined as 



B = $>kA<% At e-^* sine 
1 



kA 

xsinc 



(cj + Wfc)At 



-(ujo - w fc )At 



(53) 



Eq. (52 ) is the main result of this paper. Having a closer 



look at the definitions of A± and B, we see that pa is 
indeed time-independent, ie. Markovian. 

A relatively detailed analysis of the transition rates A_ 
and A+ can be found in App. [B] However, performing 
the integration in Eq. (B3) is only possible after choos- 



ing concrete values for the constants a.), in Eq. (21). 



We therefore restrict ourselves in the following to the 
three most prominent Hamiltonians: the rotating wave, 
the minimal coupling, and the multipolar Hamiltonians. 
Figs. [T] and [2] show A_ and A + as a function of w max - 
Since A_ and A + depend both only weakly on cv m - m , we 
assume here that uj m - ln — 0. For the minimal coupling 
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and the rotating wave Hamiltonians, A_ is essentially the Using Eq. (52 1, we find that the matrix elements of pa 



same as the usual spontaneous decay rate T in Eq. (B5|. evolve according to 



However, for the multipolar Hamiltonian, A_ is in gen- 
eral much larger than T, even for relatively small values of 
the upper cut-off frequency u; max . Fig. [2] illustrates that 
the counter-rotating terms, like cr + aj cA , in the general 
atom-field interaction Hamiltonian Hi in Eq. (27) can 
have a huge effect. For the multipolar Hamiltonian, we 
obtain transition rates A + much larger than T. Only for 
the rotating wave Hamiltonian, do we have that A + = 



(cf. Eq. (B6l) 



App. [C] shows that the master equations associated 
with the rotating wave, the minimal coupling, and the 
multipolar Hamiltonians can be written in Lindblad form 
[44j for the experimental parameters which we consider 
in this paper. This is important, since being of Lindblad 
form guarantees that the density matrices obtained from 
these equations are valid physical objects [J5J |3B] . The 
spontaneous emission of a photon can be caused by the 
de-excitation of the atomic system. But, due to the pres- 
ence of the counter-rotating terms in the system Hamil- 
tonian H\ in Eq. (20), it can also be caused by the si- 



multaneous excitation of the atom and the free radiation 
field. 



IV. FEASIBLE EXPERIMENTAL TESTS 

In this section, we calculate the stationary state pho- 
ton emission rate I ss in the absence of external driving 
for the minimal coupling, the multipolar, and the ro- 
tating wave Hamiltonians. As we shall see below, this 
rate can become enormously large, even when introduc- 
ing physically meaningful environmental response times 
At and appropriate cut-off frequencies 0J m i n and w max for 
the modes of the free radiation field. We then compare 
our predictions with the findings of actual experiments 
with single quantum dots and colour centers in diamond. 



Pn = -A + pn + A_ p 22 , 

p 12 = iw P12 + \ ( A - + A +) /°i2 + B * p 2 i , 

P21 = P*i2 ) P22 = 1 - Pll ■ 



(56) 



The stationary state of the atom (pa = 0) corresponds 
to a diagonal matrix with 



P22 = 



A_ + Aa 



(57) 



Substituting this population into Eq. (55), we obtain the 
stationary state photon emission rate 



2A^A + 
A_ +A, 



(58) 



Half of these photons have a frequency such that they 
are narrowband. This means, their frequency is close to 
the atomic transition frequency ujq. The reason for this 
is that the stationary state photon emission rate from 
the excited atomic state |2) equals A^ P22, ie. ^I S s- This 
expression is proportional to A+ (which is the rate as- 
sociated with the broadband emission of photons), since 
the off-resonant excitation of the atomic system is what 
creates the stationary state population P22 7^ in the 
first place. 

Fig. [3] shows the stationary state photon emission rate 
I ss in Eq. ( |58[ ) for the minimal coupling and the multi- 
polar Hamiltonians as a function of the cut-off frequency 
u ma i, while w m in = 0. In case of the minimal coupling 
Hamiltonian we have A + -C A_ and 7 SS is to a very good 
approximation given by 2A + . In case of the multipolar 
Hamiltonian, A + and A- are approximately of the same 
size and I ss equals A + for a wide range of cutoff frequen- 
cies. This is why Figs. [2] and [3] seem to be essentially the 
same. However, the most surprising result in Fig. [3] is 
the actual size of 7 SS . 



Detector-model dependent stationary state 
photon emission rates 



B. The minimal coupling and the multipolar 
Hamiltonians 



Having a closer look at the derivation of a general mas- 



ter equation in Section III we find that the probability 
density I(t) for the emission of a photon at time t equals 
the trace over the density matrix TZj(p ai) in Eq. (47). 



Denoting the matrix elements of the density matrix pa 
such that 



Pa 



P11 P12 
P21 P22 



wc find that 



I(t) = A_ p 22 (t) + A + p n (t) . 



(54) 



(55) 



We now have a closer look at realistic experimental 
parameters. A typical transition frequency in the optical 
regime is ujq = 3.7 • 10 15 Hz. This frequency corresponds 
to the wavelength of 500 nm. A possible estimate for the 
upper cut-off frequency for the modes of the free radia- 
tion field is w max = 3.7 ■ 10 19 Hz. This frequency corre- 
sponds to the Bohr radius [48]. As we have seen above, 
a good estimate for the environmental-response time At 
is the time it takes a photon to reach the walls of the 
laboratory. Assuming that photons travel at the speed 
of light and that the walls of the laboratory are about 
3 m away yields At = 1 • 10~ 8 s. Moreover, we assume in 
the following that T = 1 • 10 7 Hz. Using these estimates 
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and the equations for the minimal coupling Hamiltonian, 
we obtain a stationary state photon emission rate I ss of 
about 1.5 photons per second. This rate is well below 
typical detector dark count rates of about 500 photons 
per second and more or less impossible to observe exper- 
imentally. 

Substituting the same experimental parameters into 
the equations for the multipolar Hamiltonian, we obtain 
an A + of about 4 • 10 6 photons per second. This clearly 
shows that the master equation associated with the mul- 
tipolar Hamiltonian cannot be valid. In other words, the 
multipolar Hamiltonian does not identify the free radi- 
ation field seen by a photon-absorbing environment cor- 
rectly and corresponds instead to a highly unrealistic de- 
tector model. The photon emission rate predicted by the 
master equation associated with this Hamiltonian would 
have been noticed already in standard experiments with 
trapped ions. 

One way of obtaining higher stationary state pho- 
ton emission rates even in case of the minimal coupling 
Hamiltonian is to move the photon-absorbing environ- 
ment, ie. the detector, closer to the atom. The reason 
for this is that A + scales essentially as 1/At, while A_ 
remains more or less the same when At changes. For ex- 
ample, a photon-absorbing environment 10 cm away from 
the atom implies At = 3 • 10~ 10 s and the above values 
yield an I ss of about 40 photons per second. Closer detec- 
tors yields larger 7 SS . We therefore now look for atomic 
systems with a relatively large spontaneous decay rate 
r and a relatively small transition frequency coq, like the 
ones used in recent experiments with single quantum dots 
and single colour centers in diamond. These are artificial 
atoms in the sense that each of them confines a single 
electron such that its states become quantised with level 
spacings in the optical regime. 

For example, Matthiesen, Vamivakas, and Atature [5] 
studied the fluorescence from Gallium Arsenide with a 
single layer of self-assembled Indium Arsenide at a tem- 
perature of 4.2 K and measured a lifetime of T-y = 760 ps 
for an excited electronic state with a transition wave- 
length of about 950 nm. This corresponds to a decay 
rate T of about 1 • 10 9 Hz and a transition frequency 
ljo of about 2 • 10 15 Hz. Assuming moreover that At = 
3 • 10~ 10 s, we obtain a stationary state photon emission 
rate J ss of about 8400 photons per second. Measuring 
the corresponding emission rate of about 4200 narrow- 
band photons per second is experimentally feasible and 
suggests an experimental test of the validity of the master 
equation associated with the minimal coupling Hamilto- 
nian [39] • 

Similar experimental tests of the multipolar and the 
minimal coupling Hamiltonians could be performed with 
chromium-based colour centers in diamond. Recently, 
Muller et al. [3] studied the photon emission of this sys- 
tem from an excited state with a lifetime of 1 ns and a 
transition wavelength of 710 nm at a temperature of 4K. 
These parameters correspond to us^ = 2.6 • 10 15 Hz and 
r = 1 ■ 10 9 Hz. If we assume again a photon-absorbing 
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FIG. 3: Logarithmic plot of the stationary state photon emis- 
sion rate I EB for the minimal coupling and the multipolar 



Hamiltonian as a function of u> n 



while cu n 



= 0. The 



grap hs are the result of a numerical solution of Eqs. ( 58 1 and 
(B3|. 



environment which is 10 cm or less away from the colour 
centre, the condition TAt <C 1 applies and the above 
calculations hold. In this case, the master equation asso- 
ciated with the minimal coupling Hamiltonian predicts a 
stationary state photon emission rate I ss of about 2600 
photons per second. 



C. The rotating wave Hamiltonian 

Out of all the Hamiltonians considered in this paper, 
only the so-called rotating wave Hamiltonian predicts a 
stationary state photon emission rate 



J™ = 0. 



(59) 



This suggests that this Hamiltonian identifies the com- 
ponents of the atom-field system correctly and yields the 
most accurate master equation. This master equation 
is the simplest example of the general master equation 
given in Eq. ( 52 ) , 



PA 



1W [<T 3 , p A J 
) 



1 



PA 



(60) 



since all the it^'s vanish in this case (cf. Table B. As 
one can see from Fig. [I] the transition rate A- is es- 
sentially the same as the spontaneous decay rate T in 
Eq. (B5). This means, Eq. (60) is essentially the same 



as the standard Born-Markov master equation which 
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is usually obtained after applying the RWA and cer- 
tain Markovian approximations. It is Markovian and 
of Lindblad form [33] . As pointed out by Glauber [37] , 
there are hence realistic detectors which can implement 
the photon-absorbing measurements associated with this 
equation (cf. Section III B ) . 



V. CONCLUSIONS 

This paper emphasizes that there are different ways of 
quantising the states of a two-level atom inside the free 
radiation field. Depending on the choice of gauge when 
writing down the classical Lagrangian, we obtain a dif- 
ferent Hamiltonian H. On the quantum level, different 
Hamiltonians relate to each other via a unitary operator 
R (cf. Eq. ([2])). All possible Hamiltonians H are unitarily 
equivalent and possess exactly the same spectrum. The 
only difference is that each Hamiltonian refers to phys- 
ically different components as representing the "atom" 
and as representing the "field." When deriving an atomic 
master equation, we assume a photon-absorbing environ- 
ment and take the trace over the free radiation field. Each 
representation of the Hamiltonian hence results in a dif- 
ferent master equation. 

In other words, every operator R implies a different 
photon-absorbing environment and relates to a differ- 
ent detector model. The question we ask here is, which 
Hamiltonian identifies the free radiation field seen by 
a photon-absorbing environment correctly and therefore 
yields the most accurate master equation? This paper 
tries to answer this question without having to distin- 
guish between a real and a virtual photon, whose defi- 
nition depends on how the photon has been generated. 
Instead we derive master equations for a large set of uni- 
tary transformations R{ ak } (cf. Eq. (11)) which are able 
to generate the minimal coupling, the multipolar and the 
rotating wave Hamiltonians. We then propose to com- 
pare the predictions of the corresponding master equa- 
tions with actual experimental findings. 

Our derivation of atomic master equations avoids ap- 
proximations whenever possible. The only approxima- 



tions made in Section III are the assumption of an en- 



vironmental response time At and the use of second or- 
der perturbation theory with respect to an appropriately 
chosen interaction picture. The constant resetting of the 
free radiation field into its vacuum state is well justi- 
fied in the presence of an external environment which 
selects this state via the absorption of photons (einse- 
lection) |31j . This means, our model takes into account 
that photons fly away and do not return to interact again 
with their source. In the absence of an atomic source 
but the presence of a photon-absorbing environment, the 
vacuum state is the only state of the free radiation field 
which does not change in time. Moreover, second or- 
der perturbation theory is well justified as long as At 
is short compared to the characteristic time scale of the 
effective atomic evolution. This paper avoids the RWA 



and takes the counter-rotating terms in the atom-field 
interaction Hamiltonian into account. The obtained gen- 
eral master equation contains transition rates A± and B 
which depend explicitly on At as well as on the cut-off 
frequencies io m - m and w max for the modes of the radiation 
field. A good estimate for At is the light travel time to 
a detector or the walls of the laboratory. When the light 
travel time is relatively short, the time it takes a photon 
to be absorbed by the environment and the environment 
to re-thermalise should be added. 



In Section |IV[ we calculate the stationary state pho- 
ton emission rate I ss in the absence of external driving 
for the master equations associated with the minimal 
coupling and the multipolar Hamiltonian. For the ex- 
perimental parameters of recent experiments with single 
quantum dots [3] and chromium-based colour centers in 
diamond [4] , we obtain enormous emission rates (cf . Sec- 
tion IV B[ ) which would have been noticed already. Out of 
all the atom-field Hamiltonians considered in this paper, 
the only one with J ss = (cf. Eq. (59 1 ) is the so-called ro- 
tating wave Hamiltonian. Our analysis suggests that this 
Hamiltonian might yield the most accurate master equa- 
tion. The name rotating wave derives from the fact that 
this Hamiltonian is almost the same as the atom-field 
Hamiltonian obtained after applying the RWA. However, 
the rotating wave Hamiltonian is not the result of an 
approximation but the result of a specific choice of the 
unitary transformation R in Eq. Q. 

When the Hamiltonian Hi in Eq. (20) contains 



counter-rotating terms, we obtain a non-zero stationary 
state photon emission rate I ss which increases with 
decreasing At. This effect can be associated with the 
presence of a "cloud" of virtual photons localised around 
the atom. A virtual photon is a photon generated by 
a term in the Hamiltonian which cannot conserve the 
free energy, ie. a counter rotating term. When moving 
the photon-absorbing environment closer to the atom, 
the probability for the absorption of a virtual photon 
increases. Absorption of such a photon leaves the atom 
in an excited state, so that it may then decay via a (real) 
photon emission. In general, the {a/.} can be chosen 
to generate any value of J ss whatsoever. In this way, it 
becomes possible for our model to reproduce detector 
dark count rates and finite temperature effects. 
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Appendix A: Generalised PZW transformation 



In order to derive the general Hamiltonian H' 



R 



{ ak y -f/min R\ ak \, we use i n the following the Baker- 
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Campbell-Hausdorff formula 



e x Ye' x = ^2[X,Y] 



with [X,Y] r 



[x, [x, r] m _! 



closer look at H m i n in Eq. (10) and R{a k } in Eq. (11) 
shows that this calculation requires the commutators 



and [X,Y] 



(Al) 



Y. As 



[A W} (0)-r,p 2 ] 
[A {ctfc} (O)-r,A(0)-p] 

A {Qfc} (0) ■ r, t^k aL«kA 

kA 

[A w (0)t,A W} (0)- P ] 
[A {Qfc} (O).r,n {afc} (0)-r] 



2i/iA {ctfc} (0) • p, 
ifiA w (0).A(0) 



n {ctfc} (o)-r : 



ift|A {afc} (0)| 2 

m 
1? 



a 2 Irl 2 . 



(A2) 



kA 



All other commutators which occur in the calculation of 



H' equal zero. Using Eq. (A2), we obtain a Hamilto- 
nian of the same form as the Hamiltonian in Eq. ([lj. Its 
components can be found in Eq. (13). 



Appendix B: Analysis of the transition rates A± 



The aim of this section is to determine the transition 



rates A± in Eq. ( 48 ) as functions of the environmental 



response time At and of the typical cut-off frequencies 
and w max for the modes of the free radiation field. 



Using Eqs. (21) and (48) and performing the summation 



over A, one can show that 



£l<7kA| 2 
kA 



e wo \ - 



' kA 



l-|d-kl 



(Bl) 



where d and k are unit vectors in the directions of d and 
k, respectively. Here we are especially interested in the 
large volume limit, where L — > oo. Next we choose a 
coordinate system such that d points in the z-direction 
and d • k = cos ■& and replace the summation over k by 
an integration, 



kA 



L 3 
8tt 3 c 3 



dwfc w 



77 



l I M sin / dip, (B2) 
o Jo 



2tt 



where c denotes the speed of light. Taking this into ac- 
count, we then find that 



2T 



TTUJ^At 



du k f±{uJk) sin 2 



(W ± Lo k )At 



with f±(u>k) defined as 



f±(uk) 



(UJ ± LO k ) 2 



(B3) 



(B4) 



H 



/-(«*) /+(Wfc) 



H vo t 



"0 + "fe 



(u> +u; k ) 2 



TABLE I: The coefficients uf in Eq. (21 1 and the coefficients 



f±(u>k) in Eq. (B4l for the rotating wave, the minimal cou- 
pling, and the multipolar Hamiltonian. 



and with the spontaneous decay rate T defined as [50] 



(B5) 



Table [T] summarises the and f±{uJk) for these three 



cases. 



Fig. [4] illustrates a removable singularity of the coeffi- 
cients /_(u?fc) at = ujq. In case of the rotating wave 
and the minimal coupling Hamiltonian, f-(u>k) tends to 
zero when LUk tends to infinity. This means, when per- 



forming the integration in Eq. ( B3 1 , we are likely to ob- 
tain a transition rate A_ which depends only weakly on 
the cut-off frequencies w m i n and u; max as long as both are 
sufficiently different from uiq. This is confirmed by Fig. [I] 
which shows that A_ is in both cases esse ntial ly the same 
as the spontaneous decay rate T in Eq. (B5). This ap- 



plies for a wide range of environmental response times At 
and a wide range of realistic cut-off frequencies w m i n and 
w max . However, for the multipolar Hamiltonian, 
grows rapidly when LUk increases. It is therefore not sur- 
prising to see (cf. Fi g. |X[ ) that the corresponding transi- 
tion rate A_ in Eq. (B3) grows with w max . 



Out of all the Hamiltonians considered in this paper, 
the rotating wave Hamiltonian is the only one for which 



the transition rate A + vanishes, ie. 



A 



(B6) 



This applies, since / + ( w fc) = in this case, as Table 
[I] shows. The calculation of the transition rate A + for 
the minimal coupling and the multipolar Hamiltonians 
too is relatively straightforward, since f + (ujk) has no sin- 
gularity for finite values of Uk (cf. Fig. [5]). Relatively 
simple but also very accurate result s ar e obtained, when 



replacing the sin function in Eq. ( B3 ) by 75 . This ap- 



proximation is well justified, since At » 1/u)q. Using 
this approximation and the values for f + (uik) in Table [ij 
we find that 



A 



r 



nujoAt 



— h In a; 
x 



(B7) 



in case of the minimal coupling Hamiltonian, while we 
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FIG. 4: Logarithmic plot of the function /_(wfc) in Eq. (B4| FIG. 5: Logarithmic plot of the function in Eq. (B4| 



for the rotating wave, the minimal coupling, and the multi- 
polar Hamiltonian. 



for the minimal coupling and the multipolar Hamiltonian. For 
the rotating wave Hamiltonian we have /+ = 0. 



obtain 



A, 



r 



TTUjQAt 



3cc 



3 In x 



W max / UJQ + 1 



W m i n /wo + l 



(B8) 



in case of the multipolar Hamiltonian. Both rates de- 
pend strongly on the environmental response time At as 
well as on the cut-off frequencies w m ; n and w max . This 
is illustrated in Fig. [2] which shows A + as a function of 
w max for a fixed At and w m ; n = 0. 



the form 

PA 



-1W [CT3,PA 
2 



E^ 



UPkL\ - -{L\Li,p A } 



(C2) 



can now be achieved via a diagonalisation of M. To do 
so, we write M as 



M 



£>|Ai)(A<l> 



(C3) 



Appendix C: Lindblad form of the general master 
equation 

In this appendix we diagonalise the general master 
equation in Eq. ( 52 ) and identify necessary conditions 
for this master equation to be of Lindblad form ffl] . To 



do so, we notice that Eq. (52 1 is an equation in first stan- 
dard form 



Pa = -iw F3,AaJ 

2 



n,m— 1 



VnPAvL ~ -^{^m^n, Pa) 



, (CI) 



where {, } denotes the anticommutator and where the 
M nrn are the matrix elements of a Hermitian operator 
M. A transfer of this equation into Lindblad form, ie. in 



where the and the |A$) are the eigenvalues and eigen- 
vectors of M. One can now easily check that this matrix 
is diagonalised by the unitary matrix 



U 



Ei e *>< A * 



(C4) 



where the |e<) are the canonical vectors |ej) = (1,0) T 
and fa) = (0, 1) T . The matrix U is indeed such that 

&MU = ( ^ A ° 2 ) • (C5) 
Now, defining new operators Li such that 



CTl 
<T 2 



(C6) 



and substituting them into Eq. (CI) gives the diagonal 



master equation of Lindblad form in Eq. ( C2 ) 
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is therefore of Lindblad form with one of the two Lind- 
blad operators being zero. As Eq. (60 1 illustrates, the 



remaining non-zero Lindblad operator is the bare atomic 
lowering operator er~. 



To check whether Eq. ( C8 ) applies to the master equa- 



tions associated with the multipolar and the minimal 
coupling Hamiltonian, we now proceed as in App. |B"]and 
calculate B as a function of At, ujq, w max , and L. To do 
so, we write the normalised atomic dipole moment d and 
the general wave vector k as 




and k 



sin ■& cos tp 
sin $ sin tp 
cos^ 



(C9) 



where i? and tp are the usual polar coordinates for real 
vectors in three dimensions. Two normalised polarisation 
vectors eti and are then given by 



FIG. 6: Logarithmic plot of an upper bound for A+ A_/|i3| 2 
for the multipolar and the minimal coupling Hamiltonian ob- 
tained from a numerical integration of Eqs. (B3l and (531 
after taking into account that \di + d| + 4\ 2 < 1 and that 
cos 2 (o;oAt) < 1, while assuming Uxnin = 0. The figure con- 
firms that Eq. ( |C8[ ) holds for the experimental parameters 
which we consider here and that the corresponding master 
equations are in general of Lindblad form. 



COS # cos tp 
eki =| cos $ sin tp | e k2 = 
— sini? 




(CIO) 



One can easily check that these two vectors are pairwise 
orthogonal and orthogonal to k. Substituting these into 
Eq. ( 53 1 and performing the spatial integration over & 



and tp, while taking the definition of the coupling con- 
stants gkA in Eq. (21 1 into account, we find that 



We now turn our attention again to the general master 
equation in Eq. (52 1. This master equation is of the same 
form as Eq. 
<r + and (72 - 



(CI ), if we define the operators cr, by a\ 
o~ , while defining the matrix M by 



M = 



A+ B 
B* A_ 



(C7) 



The corresponding Lindblad operators Li are always 
traceless. This means, the corresponding master equa- 
tion \G2\ is of Lindblad form, if the eigenvalues \i are 
either positive or equal to zero. This applies when 
det(M) > 0, ie. when 



A>A_ > \B\ 



(C8) 



The definitions of the coefficients A± and B can be found 
in Eqs. (|48l and p3|. 



As it has been pointed out already in App.|Bj to calcu- 
late A± and B a concrete representation of the atom-field 
Hamiltonian H must be chosen. In the following, we pro- 
ceed therefore as above and restrict ourselves again to the 
rotating wave, the minimal coupling, and the multipolar 
Hamiltonian. For example, in the case of the rotating 
wave Hamiltonian, we have ul = (cf. Tab. |l| which 
implies A_ = B = 0. This means, one of the eigenval- 
ues of M is zero and the other one is positive. Eq. ( 52 ) 



B = — 



2L 



ttujq At 



-iujq At 



(4 + 4 + df) 



dujj 



u i u k 

2 2 



sin 



(lu q + u> k )At 



-{ojq - uj k )At 



(Gil) 



Using trigonometric relations and approximating the in- 
tegral over cos(w£,A£) by zero, this expression simplifies 
to 



B = 



At 



cos(w At) [d\ 



4 



4) 



dujk 



u\ u t u k 



(C12) 



Concrete expressions for the can be found in Table [i] 
Replacing 1 4 + 4 + 4 1 an( l I cos(wo At) | by their respec- 
tive maximum of one, we obtain an upper bound for \B\. 
Fig. [6] uses this bound and verifies Eq. ( C8 ) for the mul- 
tipolar and the minimal coupling Hamiltonian for a wide 
set of experimental parameters. This shows that the mas- 
ter equations associated with these two Hamiltonian are 
in general of Lindblad form. 
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